Generalized Liischer Formula in Multi-channel Baryon-Meson Scattering 
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Liischer's formula relates the elastic scattering phase shifts to the two-particle energy levels in a 
finite cubic box. The original formula was obtained for elastic scattering of two massive spinless 
particles in the center of mass frame. In this paper, we consider the case for the scattering of a 
spin 1/2 particle with a spinless particle in multi-channel scattering. A generalized relation between 
the energy of two particle system and the scattering matrix elements is established. We first obtain 
this relation using quantum-mechanics in both center-of-mass frame and in a general moving frame. 
The result is then generalized to quantum field theory using methods outlined in Ref. [1]. We verify 
that the results obtained using both methods are equivalent up to terms that are exponentially 
suppressed in the box size. 



I. INTRODUCTION 

Low-energy hadron-hadron scattering plays an impor- 
tant role for the understanding of strong interaction. 
However, due to its non-perturbative nature, it should 
be studied using a non-perturbative method like lattice 
Chromodynamics (lattice QCD). Lattice QCD can tackle 
the problem from first principles of QCD using numer- 
ical simulations. By measuring appropriate correlation 
functions, energy eigenvalues of two-particle states in a 
finite box can be obtained. Liischer found out a relation, 
now commonly known as Liischer's formula, which relates 
the energy of two-particle state in a finite box of size L, 
E(L), to the elastic scattering phase 8(E(L)) of the two 
particles in the continuum [2-5] . While the former could 
be obtained in lattice QCD simulations, the latter fully 
characterizes the scattering property of the two particles 
and could in principle be measured in corresponding ex- 
periments. Thus, this relation opens up the possibility of 
lattice study of hadron-hadron scattering. 

Liischer's formalism has been utilized in a number 
of lattice applications, e.g. linear sigma model in 
broken phase [6], hadron-hadron scattering both with 
quenched approximation and unquenched configuration 
that contains dynamic quarks [7- 17]. However, the origi- 
nal Liischer's method is restricted to elastic scattering of 
massive, spinless particles in the center-of-mass (COM) 
frame of the two particles. Some of these constraints 
restrain the applicability of the formalism to general 
hadron scattering. For example, in real simulations, one 
is usually restricted to only a few lattice volumes and 
since the energies in the box is quantized, and taking into 
account the fact that excited energy states are more dif- 
ficult to measure numerically, one ends up with a rather 
poor energy resolution when compared with the exper- 
iments. To overcome this difficulty, one can of course 
consider asymmetric volumes [18-20], or boosting the 
system to a frame that is different from COM [21-25], 
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both will enhance the energy resolution of the problem. 
Another possible generalization is to use the so-called 
twisted boundary conditions advocated in Refs. [26-29]. 

Generalizations to particles with spin [30-32] is also 
possible. For example, in Ref. [33, 34], Liischer's formula 
has been extended to elastic scattering of baryons. 

It is also important to extend Liischer's formula to the 
case of inelastic scattering which is commonly encoun- 
tered in hadronic physics. Attempts have been made 
over the years, see Refs. [1, 35-41]. 

In this paper, we would like to synthesize the above 
mentioned generalizations by trying to seek a formula 
(or a formalism) that is applicable for multi-channel scat- 
tering of particles with spin in a possibly moving frame 
(MF). For this purpose, we consider two-particle to two- 
particle scattering processes in which one particle has 
spin 1/2 (we will henceforth call it a "baryon") while the 
other particle remains spinless (we will call it a "meson" ) . 
The scattering process we study could be a multi-channel 
scattering beyond a certain threshold. But in each chan- 
nel, the two-particle states (initial or final) always has 
the feature that one particle has spin 1/2 while the other 
is spinless. We will comment on possible extensions to 
this restriction in Sec. IV. As it turns out, the basic for- 
mulae obtained are the same in non-relativistic quantum 
mechanics and in massive quantum field theories, we will 
start our discussion in the former case. Generalization to 
quantum field theory can be achieved by using methods 
outlined in Refs. [1, 34] which is detailed in Sec. III. 

The organization of the paper is as follows: In Sec. II, 
we start out by discussing a quantum-mechanical model. 
Liischer's formulae are obtained first in the case of COM 
frame in subsection II A and then generalize to moving 
frames in subsection II B. In Sec. Ill, we generalize the 
results obtained in Sec. II to quantum field theory. This 
is achieved first in the single channel situation and then 
to the two-channel scenario. We then compare the results 
obtained in quantum field theory with those from quan- 
tum mechanics in Sec. II. It is shown that they are in fact 
equivalent up to terms that are exponentially suppressed 
in the large volume limit. In Sec. IV, before we conclude, 
we will also discuss possible applications of our formulae 
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in real lattice simulations and comment on some possible 
extensions in the future. 

Some calculation details arc summarized in the appen- 
dices. To be more specific, single channel scattering for 
particles with spin in infinite volume are reviewed in ap- 
pendix A. For reference, single channel Liischer's formu- 
lae are also provided in this appendix. Calculations of 
loop summation/integration in the case of quantum field 
theory are summarized in appendix B. 



II. LUSCHER'S FORMULA FOR 
BARYON-MESON SCATTERING IN 
NON-RELATIVISTIC QUANTUM MECHANICS 
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Note that due to spin degrees of freedom, each compo- 
nent is still a two-component spinor. At large r where the 
potential V(r) vanishes, the wave functions of the scat- 
tering states can be chosen to have the following forms: 
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Liischer's formula in COM frame 



1. Two-channel scattering in the continuum 
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In this section, using non-rclativistic quantum mechan- 
ics in COM frame, we will briefly discuss two-channel 
potential scattering in the continuum for the case of two 
stable particles: one with spin (a meson) and the other 
with spin i (a baryon). We will follow the discussion 
in Rcf. [35]. The potential V(r) between the two parti- 
cles is assumed to have finite range, i.e., V(r) =0 with 
r = |r| > R for some positive R, but the potential itself 
could in principle be spin-dependent so that the spin of 
the fermion might change during the scattering process. 
We assume that there exists a threshold Et > and the 
energy of the two-particle system becomes 



2^2 
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where \i\ and [12 are the reduced mass of the two-particle 
system below and above the threshold, respectively. In 
the COM frame, one only has to denote the momentum of 
one of the two particles. For definiteness, we denote the 
momentum of the fermion as ki and 1<2 in the first and 
second channel, respectively. The magnitude of them are 
ki = |ki| and &2 = |lt2 J • Obviously, for energies below 
the threshold &2 will become pure imaginary. 

The wave function of two-particle system, after fac- 
toring out the trivial COM coordinates, has a two- 
component form in the case of two-channel scatter- 
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In the above expressions, Xs is an eigenstate of spin 
angular momentum s z of the baryon with eigenvalue 
s = -1/2,1/2. The scattering amplitudes like M^Jf] 

depends only on the corresponding angles ki f and k2 f . 
In order not to confuse with the matrix M introduce in 
quantum field theory in Sec. Ill, we have added a su- 
perscript (NR) to stand for the case of non-relativistic 
quantum mechanics. The notation in this paper is as fol- 
lows. Subscripts like i and j refer to the channel, and 
take the values 1 or 2. It is seen that, in the remote past, 
■0i;s( r ) becomes a pure incident plane wave in the first 
channel with definite linear momentum ki and definite 
spin s. Similarly, in the remote past, ■02 ; s(r) represents 
an incident plane wave in the second channel with definite 
momentum k2 and definite spin s. It is also clear that 
the above wave functions, ipi- tS (r) with s = —1/2, 1/2 are 
linearly independent and in fact they are also complete 
in the sense that any eigenfunction of the Hamiltonian 
must be a linear superposition of them. 

As is said, the scattering amplitudes M^^l for s',s = 
1/2,-1/2 depend only on the corresponding angles and 
can be expanded into spherical harmonics, see for exam- 
ple Ref. [42]. For simplicity, we have chosen the z-axis 
to coincide with the incident momentum, ki and k2 for 
i/>i ;s (r) and ^^(r), respectively. The scattering ampli- 
tudes then takes the following form: 



< s ^(k 2 .r) 



1=0 j=i-h 



J2 E VM2l + l)SfH% ms ,S l } M . 0s Y lm (r) . 
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In the above expressions, S l / M . ms , = (JM\lm; ^s') and Sja/-o s = (JM\10; |s) are the Clebsch-Gordan (CG) co- 
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efficients. 

The quantities Sfj" are the .S-matrix elements where J 
is the quantum number of the total angular momentum. 
Since we are dealing with the scattering of a spin 1/2 
particle, J can take two possible values for a given I, 
namely J = Z± 1/2. To simplify the notation we will also 
denote the corresponding 5-matrix elements by 



S 



J=l±l/2,l 
ij 



s 



l± 



(7) 



Needless to say, one also has similar expressions for 

M^.^l and M^.^l- I n the following, we will also need 
the so-called spin spherical harmonics defined as 



2. Two-channel scattering in a cubic box 



Now we put the two-particle system into a cubic box of 
size L and impose the periodic boundary condition. The 
potential then becomes Vt(r) = J2 n + n -^)' n £ 
We divide the whole space into two regions: the inner 
region and the outer region. In the inner region, every 
point satisfy the condition: |r+nL| < R, for some n e T? 
while in the outer region, ft — {r| : |r+n£| > R, n € Z 3 }. 
In the inner region, the solution to Schrodingcr equation 
of the system is 



JM:ms ' 



(8) 



which is an eigenfunction of the total angular momen- 
tum of the system. The above expressions are just di- 
rect generalizations to the spin-dependent single-channel 
scattering. For convenience, some relevant formulae are 
collected in appendix A. 

With the spin spherical harmonics defined in Eq. (8), 
we can expand the wave function in the following form: 

fe(r) = ]T \Z^W+rjS% 0s Wi^r)Y^ M (r) , (9) 

JMl 

where the radial wave functions of Schrodingcr equation 
are denoted by W^ji(r). In the large r region, they have 
the following asymptotic forms 



Wx.,ji{r) = 



2ik\r 



l+l e -ikir] 
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W 2 , Jl (r) = 



2fr\ 



I Ml oJl^.ik-1 r 



Jl„ik 2 r 
22 e 



+ (-1) 
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It is obvious that two radial wave functions W\-ji(r) 
and W2-ji(r) are linearly independent. Since the radial 
Schrodingcr equation has two linearly independent solu- 
tions which are regular at the origin, denoted by Ui-ji(r), 
these radial wave functions can be expressed as linear 
superpositions of two radial wave functions Wi ; j/(r) [35]. 



*w = £ E E 

1=0 J=l-l M=-J 



1 ^2,h;JMlUj-Jl{r) 
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Since Ui-ji(r) can be expressed as linear superpositions 
of two radial wave functions Wi-ji(r), in the outer region 
the wave function is given by 



V>(r) 



OO 

E E 

1=0 J=l- 



J 

E 

I M=-.J Li 
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with Fi-jMi being some non-trivial constants. On the 
other hand, in the outer region f2, the wave function 
V>(r) is also a linear superposition of the so-called sin- 
gular periodic solutions [4] to the Hclmholtz equation, 
Gi ; jj\/i(r; kf), also known as the Green's functions. Thus 
we also have 



( oo l+k J \ 
E E E Vl;JMlG Vl jMl(r;kf) 
l=0j=l-^M=-J 

oo l + k J 

E E E V2-jMiG 2 ;.jMi{r;kl) 

\ l=0j=l-iM=-J J 



(14) 



From Rcf. [4, 34], the Green's function for particles with 
spin takes the following form: 



G 4; jMz(r,/c l 2 ) 



(-i)'fc; 

47T 
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l'=0 J'=l 



E 



j' 
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I M' = -,V 



M 

i;JMl\J'M'V 



(15) 



where the explicit form of M.i-jMl;J'M'V (kf) can be found tutc k 2 — kf into Eq. (All) to get M^jmi-.j' M'i'(k 2 ). 
in appendix A. For M.i-jMl;J'Wl'(kf), one has to substi- 
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To simplify the notation, wc define F\.jmi = 



and V2;jMi = \fJF\V2;jMi- Then the equivalence of 
Eqs. (13), (14) leads to four linear equations involving 
F- j mi's and Vi-jMi's, 



F\-jmi{S'1[ + 1) + F 2 -jmi 
-iFi ; jMi(Sii — 1) — iF 2 -jMi\ 

F2;JMl{S'22 + 1) + Fi-JMI" 

-iF 2 -jMi(S22 — 1) ~ iF\-jMi^ 



fc l oJl 
D 12 



= E m 
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fc 2 a Jl 



«2 Q Jl 
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J'M'l' 



47T 



-V2;J'M'V 



(-l)'fc 
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4- 



V2-JMI 



One can eliminate the coefficients v\-jmi and V2-,jmi cas- Using matrix -M^jmvj'M'V' we can wrl te the secular 
ily, leaving behind a set of linear equations for Fi-jmi equation for Eq. (16) as 
and F 2 -jmi- In order to have non-trivial solutions for 
them, the determinant of the corresponding matrix must 
vanish. Let us now define 

^i-JMhJ'M'l'( k i) = ■M-i;JMl;J'M'l'{k 2 l ) ±i8jj,8 M M'5u> . 

(17) 



\;JMl;J'M'V 
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/V2 
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We will also use the more compact matrix notation 
M.12 — Ali 2 i i- Assuming that the matrix Al^ is 



non-singular, we define a unitary matrice f/, as 

, -l" 



i-JMhJ'M'V — 



Mf {M 

f M l + i 
\Mi-i 



JMhJ'M'V 



(19) 



JMhJ'M'V 



With this matrix, the generalized Liischer's formula for 
two-channel baryon-meson scattering may be written in 
an equivalent form 



Ui-JMl-J'M'V — SuSjj'SmM'Sw 
^S&jj'SmI 



^S^Sjj'Smm'Sw 



2;JMhJ'M'V 



S^JJ'SmM'&U' 



(20) 



3. Liischer's formula with definite cubic symmetry 

In the sector with definite cubic symmetry, the basis oi 
the representations are |r, £; J, l,n), where T labels the 



symmetry sector (the irreducible representation of the 
cubic group), £ runs from 1 to the number of dimension 
for the irreducible representation, and n runs from 1 to 
the multiplicity of the irreducible representation. Then 
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it can be expressed by the linear combination of \JMl), 
where the matrix Ai is diagonal with respect to T and £ 
by Schur's lemma. In the symmetry sector T, the general 
formula (20) is reduced into: 



I k 2 cJl 



M U 2 (T) 



qJl 

"19 



(21) 



where Ui(T) represent a linear transformations of the 
space Ha(T) with J < A. In terms of matrix element, as- 
suming that the irreps T appears only once, we will label 
them as Ui(T) jh'-ji- To write out a more explicit for- 



mula, we should consider the definite cubic symmetries. 
For the case of half-integer total momentum J, we need 
to consider the double cover group of O denoted by O d , 
which contains 48 elements and can be divided into 8 
conjugate classes: A\, A 2 , E, Ti, T 2 , Gi, G%, H. For in- 
stance, for J = |, |, |, |, A = 4, the decomposition into 
irreducible representation are given by \ = G\, § = H, 
|=i70G 2 ,|=i?0Gi0G 2 respectively [43]. 

Now we focus on the G\ and G 2 sector. In Gi sector, 
there is a mixing between J = i (/ = 0, 1) and J = |(£ = 
3, 4). If we neglect this mixing, then there is only mixing 
within J = 1/2 between I = and 1 = 1, i.e. between s 
wave and p wave. In this case, Liischer's formula takes 
the following form 





U 



Do 



= 



(22) 



In G 2 sector,the situation is similar, and there exists 
a mixing between d-wave (/ = 2) and /-wave (I = 3). 
Liischer's formula takes exactly the same form as Eq. (22) 
except that all the labels of I', I = 0,1 are replaced by 
l',l = 2,3 and ,7=1/2 by J = 5/2. 

B. Liischer's formula in moving frames 

In this subsection, we extend two-channel Liischer's 
formula that has been obtained in the previous subsec- 
tion for meson-baryon scattering to moving frames (MF). 
This is necessary in some lattice applications since it pro- 
vides more low-momentum modes for a given lattice. Al- 
though we will only focus on the case of meson-baryon 
scattering, similar steps can be followed in the case of 
hadron scattering with arbitrary spin. We will follow the 
notations in Ref. [21] below. 

We denote the four momenta of the two particles in the 
lab frame, which is the frame in which periodic boundary 
conditions are applied, by 

k = (E 1 ,'k), P-k = (£ 2 ,P-k) , (23) 

with E\ = ^/k 2 + m\ and E 2 = \J (V — k) 2 + being 
the energies of the two particles in the lab frame and mi 
and m 2 being the mass value of the baryon and meson 
respectively. The total three momentum P / of the 
two-particle system is quantized by the condition P = 
(27r/L)d with d G Z 3 . The COM frame is then moving 
relative to the lab frame with a velocity 

v = P/{E 1 +E 2 ) . (24) 



In the COM frame, the momenta of the two particles will 
be denoted by k* and (— k*), respectively, k* is related 
to k by conventional Lorentz boost: 

k*"= 7 (kll-v^ 1 ),k* ± =k ± , (25) 

where the symbol _L and || designates the components of 
the corresponding vector perpendicular and parallel to 
v, respectively. For simplicity, the above relation is also 
denoted by the shorthand notation: k* = 7k. A similar 
transformation relation holds for the other particle. 

Let (j)(r) represent the wave function of the system in 
lab frame where r is the relative coordinate between the 
two particles. Next, we enclose the system in a cubic 
box with finite size L > 2R and apply periodic boundary 
conditions to 0(r). On the other hand, this wave function 
can be related to the COM wave function <p CM (v) by a 
Lorentz transform. Periodic boundary conditions in 4>(t) 
then implies that (j) (r) fulfills the so-called d-periodic 
boundary condition [21, 23, 24]: 

(j) CM (r) = e" aA - a (j) CM {r + 711L) (26) 

withn G TI? , a = l + (m\-ml) / 'E* 2 , E* = \Jm\ + k* 2 + 
\J m 2 , + k* 2 is the total energy in the center of mass frame 
and d = (Z/27r)P. 

In two-channel scattering, the COM wave function of 
the system can be written as: 

where the form of 0f M (r) is the same as ^(r) given in 
Eq. (2) in subsection II A. In the outer region, </)f M (r) 
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can be also expanded in terms of modified Green's func- tion Gf. JMl (r; k* 2 ) 



J 



0° M (r) 



EE E Vi;JmGl. Jm (T;kf) X 

!=0j=i-iM=-J 

EE E ^;jmGi. JMl (v;kf) 



(28) 



I 

Just like in COM frame, the corresponding Green's func- tions Gf. JMl (r, k* 2 ) are given by an analogous expansion, 
I 



(-l)'fc 



4~ 



OO ^ +^ 

E E 

l'=0j>=l>- 



J' 

E 

I M' = -.T 



M 



i;.JMl:J'M' 



, r (fc* 2 )y;4,(f)jHfc*r) 



(29) 



where the explicit expression for Aif. 
found in appendix A, e.g. 
unitary matrix as 



i-.JMhJ'M'V 



can be 



Eq. (A19). We also define a 



+ i 



i-.JMUJ'M'V 



\Mf- 



(30) 



JMUJ'M'V 



Then, following similar steps as in previous subsection, 
we can also obtain Liischer's formula in MF which takes 
exactly the same form as Eq. (20) except that all the 
matrix elements of Ui are replaced by those of Iff. 

In moving frames, to describe the scattering phase in 
definite symmetry sector, one should consider the cubic 
lattice space group Q which is the semi-direct product of 
lattice translation group T and the double-covered cu- 
bic group D . The representation of Q can be char- 

I 
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acterized by two indices: total three-momenta P and a 
representation T of the little group corresponding to mo- 
mentum P. As examples, we will only discuss moving 
frames with total momentum P = (^)e^ (MF1) and 
with P = (^)(ei + e 2 ) (MF2) [22] in the following. To- 
gether with MF1 and MF2, another moving frame with 
P = (27r/L)(ei + e 2 + e$) (MF3) has been discussed in 
Rcf. [25] for the case of single- channel (elastic) scatter- 
ing. Strategies for the construction of lattice operators 
are also analyzed. In principle, these results could be 
generalized to the case of multi-channels as well. 

In the case of MF1, the little group is C% v , which has 
7 conjugate classes: Ay, A2, B-y, B2, Ey, E2, £3. In Ey 
section, there is mixing between J ~ \ and J = \ [44] . If 
we neglect this mixing and assume an angular momentum 
cutoff of A = 4, Liischer's formula becomes 





U 2;i 0;i 




U 2;i 0;i,l 



U. 



24,1,4,1 



-S. 



22 



(31) 



In the case of MF2, the little group is Cy v , which has 
5 conjugate classes: Ay, A2, By, B2, E. When A = 4, in 
E section, there is mixing among J = ^,J= | , J = | 
and J = I [44] . The formula becomes quite complicated. 
However, if we only consider the case of J = |, the for- 
mula has the same form as Eq. (31). 



III. GENERALIZED LUSCHER'S FORMULA IN 
QUANTUM FIELD THEORY 



In this section, we will describe the generalized 
Liischer's formulae for meson-baryon scattering in quan- 
tum field theory We will follow the strategy outlined in 
Refs. [1, 45], sec also Ref. [46, 47]. In what follows, we 
will first perform the discussion in a single channel case. 



It is then generalized to two-channel case in a straight- 
forward manner. Then we compare what we obtained in 
quantum field theory with the results obtained in non- 
relativistic quantum mechanics in the previous section 
and show that they are equivalent, apart from possible 
corrections that are exponentially suppressed in the large 
volume limit. 

The discussion using relativistic quantum field theory 
has an advantage, namely the results are easily trans- 
formed to any frame, the COM frame or the lab frame 
(moving frame), both of which have been discussed in 
the previous section. In the single channel scenario, we 
denote the masses of the meson and the baryon as mi 
and rri2, respectively. The total four momentum of the 
two-particle system is denoted as P — (E, P) and the 
quantities in the COM frame will be denoted by adding 
a * to the corresponding quantity. Thus, for example, the 
COM frame four momentum is denoted as P* = (E* , 0). 
The individual three-momentum of the two particles will 
be denoted by q* and q*, the magnitude of which 
(q* 2 = q* 2 ) being: 



4q* 2 = E* 2 



2{m{+m 2 2 ) + 



{ml 



E 



*2 



(32) 



A. Single channel case 

We start by deriving an expression from quantum field 
theory in the single channel case. We first use the method 
that have been studied in [1, 45] to obtain the quanti- 
zation condition, based on the periodic boundary con- 
ditions hence the total momentum being P = (2ir/L)n 
with (n G Z 3 ). 

The basic idea in Rcf. [1, 45] is the following: two- 
particle spectrum of the system in a finite box can be 
determined from the poles of an appropriate correlation 
function in the energy plane. Thus, one defines 



C(P) = 



J(-P-x+Ex" 



\a(x)aH0)\0) , (33) 



where P = (E, P) is the total four-momentum of the 
two-particle system. The generic interpolating operator 
cr(x) is chosen to have an overlap with the two-particle 
states that we are interested in (in our case, a meson and 
a baryon) and J L . = J L d x stands for the space-time in- 
tegration over the finite volume. Two-particle spectrum 
are exactly those poles in the E plane of C(P). 

The correlation function C (P) may be built from a se- 
ries of contributions illustrated by diagrams in FIG. 1(a). 
In this figure, a solid line with a black dot stands for a 
full propagator of the baryon while a dashed line with a 
dot stands for that of the meson. A circle with a symbol 
iK inside represents the Bcthc-Salpeter kernel which con- 
sists of all amputated two-particle irreducible diagrams. 
A circle with a symbol a or denotes the interpolating 
operators in Eq. (33). Using iK we may write the C(P) 
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FIG. 1. Horizontal dashed lines represent the meson prop- 
agators while the solid lines represent the baryon propaga- 
tors. The corresponding full propagators are represented by 
dashed/solid lines with a black dot. The dashed rectangles in- 
dicate finite volume momentum sum/integrals. Each dashed 
vertical line indicates an insertion of the factor F entering 
Eq. (43) which carries the volume dependence of interest, (a) 
The diagrams that build up the correlation function C(P) de- 
fined in Eq. (33); (b) and (c) Diagrams that build up matrix 
element A and the scattering amplitude iM, respectively; (d) 
The resulting series for the subtracted correlator defined by 
Eq. (42). 



in the following form: 

C{P)= f a q \Z x A x Z 2 A 2 \ q a\ 

+ f a q [ZiAiZ 2 A 2 ], iK M , [Z 1 A 1 Z 2 A 2 } q , a\, 

J L.a.a' 



(34) 



where we have adopted the shorthand notation for the 
two-particle propagators 

[Z 1 A 1 Z 2 A 2 ] q = [Zi(g)Ai(?)] • [Z 2 (q)A 2 (P - q)] . (35) 

Denoting the interpolating field for the baryon and the 
meson by 4>i(x) and 4> 2 {x), respectively, the full propa- 
gators appearing in the above equations read 

Z t (q)A t (q) = J d^e^iM^MO)) > (36) 
*(?' 1 7m + m i) 



Ai(<?) 
A a («) 



q 2 — m\ + is 



q 2 — m\ + ie 



(37) 
(38) 



The factors Z\ (q) and Z 2 (q) are the corresponding dress- 
ing functions for the baryon and the meson, respec- 
tively. In Eq. (34), two- particle intermediate states are 
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summed/integrated in a manner that is appropriate for 
the finite volume, namely, 



L,q 



2 71 



(39) 



The kernel iK is related to the Bcthe-Salpeter kernel BS 
as 



iK, 



1,1' 



iBS(q,P-q,-j,-P + j) , 



(40) 



with BS the sum of all amputated two-particle scattering 
diagrams which are two-particle-irreducible. Finally the 
factors a q and ctJ denote the coupling of the interpolating 
operator to the two-particle states. 

Normally the interpolating operator is the product of 
two interpolating operators for the two hadrons being 
considered, e.g. cr(x) = 0i(x)</>2(x). In order to have 
an overlap with the desired states, interpolating opera- 
tors are usually designed to carry definite J PC quantum 
numbers. Therefore, without loss of generality, we as- 
sume that a(x) carries definite parity. In our case, the 
operator cr(x) must carry Dirac indices as (j>i{x) and since 
7o is the Dirac matrix responsible for parity transforma- 
tion for Dirac spinors, we assume that cr(x) commute 
with 7q. In other words, in what follows, 70 may be 
considered as a c-number rather than a matrix in Dirac 
space. 

The kernel iK and the propagator dressing functions 
^1,2(9) have only exponentially suppressed dependence 
on the box size L. Such dependence is always assumed 
to be small and negligible in the large volume limit. The 
dominant power-law volume dependence enters through 
the discrete momentum sums in the two-particle loops 
in Eq. (34) , the details of these calculations are outlined 
in appendix B, see Eq. (B12). Basically, each summa- 
tion/integration can be split into two parts: 



I — loc 1 + IfV > 



(41) 



where loo designates the infinite volume result of the loop 
integral and Ipv contains the finite volume corrections. 
For each of the loop summations in FIG. 1(a), if we take 
loo in each loop, we then recover the infinite volume cor- 
relator C°°(P). As shown in appendix B, C°°(P) does 
not contain the two-particle poles that we are looking for. 
The part of interest is the finite volume correction: 



C FV (P) = C(P)-C°°(P) 



(42) 



where C (P) and C°°(P) are the finite and infinite 
volume correlation functions, respectively. Diagrammat- 
ically, C FV (P) is obtained by keeping at least one inser- 
tion of Ipy in the two- particle loops in Eq. (34). Let 
F be the kinematic factor associated with the factor of 
Ifv ■ These contributions are shown in Fig. 1(d), leading 



to the following general result: 



C 



FV 



AFA' + AF(iM)FA' 
1 



—AF 



-A' 



(43) 



1 + iMF 

The finite volume correction C FV contains the two- 
particle poles we are looking for. 

The factor F is discussed in some detail in appendix B, 
see Eq. (B14). It arises from the integration/summation 
of the intermediate states represented by the loop di- 
agram, see Eq. (B6). Due to the fermion propagator, 
the factor F is in principle a matrix in Dirac space, as 
shown in Eq. (B7). However, it only contains the matrix 
70, which is the matrix responsible for the transforma- 
tion of a Dirac spinor under parity. Since in practical 
applications we normally choose the interpolating opera- 
tors cr{x) to carry definite parity quantum number, this 
means that 70 can be replaced by its eigenvalues: ±1. 
In other words, we can simply treat the function C(k) 
defined in Eq. (B7) as c-numbers. Thus, in the COM 
frame, we simply define: 



F = C(q*)F,M = C(q* 
with this, Eq. (43) becomes 

C FV = -AFA 1 + AF(iM)FA' 4 



(44) 



-AC(q*)FA' + AC(q* )F(iM)FA' + 
1 



-AC(q*)F 



^A! 



1 + IMF 



(45) 



The factors A and A! appearing in Eq. (43) may be 
expressed by appropriate matrix elements of the interpo- 
lating operator, as shown in Fig. 1(b). In COM frame 
the amplitudes A and A' now read, 




(0|(r(0)|k%-k*;in)| k .| =9 
(k*,-k*;out|o-t(0)|0) 



(46) 



|k*|=g- 



Both of these two amplitudes can be viewed as two- 
component spinor in spin space. It is more useful to 
express the abstract formula (45) in angular momentum 
basis. For this purpose, we expand the matrix elements 
in terms of spin spherical harmonics defined in Eq. (8): 




■AjMiYjjjik*) 

■a'jmiYjU**) 



(47) 



where summation over repeated indices (i.e. J, M and 
I) are understood. 

Similarly, working in states with definite parity, F and 
M are viewed as 2 x 2 matrices in spin space, which can 
be expanded as well, 
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F(k*,k*') = ——FjMl; 

47T 



J'M'l 



Y 



(48) 



M(k*,k*') =47rM. 



JMl-J'M'l 



(49) 



(1 - ('U - 

Here the notation F 7 | f (k*)Y J ,l I , (k*') stands for the di- 
rect product of two spin spherical harmonics. Then the 
abstract formula Eq.(45) remains valid except that all 

I 



r 



implicit indices are in angular momentum space. For ex- 
ample: 



(AC (q* )FMFA')j 1 



= A 



M 1 h;J 2 M 2 hM J 2 Mzlr,J 3 M 3 l 3 r J 3 M 3 h-JiMil 



Ft 



MA, 



(50) 



where the matrix F 



JMl-J'M'l' 



is given by: 



F, 



JMl-J'M'V 



8irE' 



■($jj'$mm'Sw + iF t 



FV 

JMl-J'M'l 



) (51) 



with Fjmi-jim'V gi yen m Eq. (B16) in the appendix. 

The matrix M in Eq. (45) is the scattering amplitude 
illustrated in FIG. 1(b) and FIG. 1(c) which can be re- 
lated to the non-relativistic quantum mechanical scatter- 
ing matrix M^ NR > . Following the discussion as in Rcf. [5] , 
the relation is found to be 



M 



JMl-J'M'V 



= 8wE*M 



(NR.) 

JMl-J'M'l' ■ 



(52) 



where M 



{NR.) 

JMl-J'M'V 



given by Eq. (A4) in appendix A. 



The quantization condition can now be obtained for 
C FV (P) which manifest itself as a series of two-particle 
poles. This means that the matrix between A and A' 
must have divergent eigenvalues hence satisfy the follow- 
ing condition: 



dct(l + iMF) = 



(53) 



This is the so-called quantization condition for the two- 
particle poles in a finite box. In order to compare this 
with the conventional Liischcr's formula obtained by non- 
relativistic approach, wc substitute the definitions of 

I 



I 

Mjmi-j'M'I' and Fjmuj'M'I' in Eq. (52) and Eq. (51) 
into Eq. (53). We then obtain the following result: 



Aet(tan 5 jiF j j, Mn , - 5jj>6mm'Sw) = • (54) 



Using an equality which will be shown in appendix B, 
namely Eq. (B19), it is seen that F is related to the 
matrix A4 d up to a possible phase. In the case of strong 
interaction, the phase factor does not enter, since the 
Hamiltonian of the system conserves parity and there 
can be no scattering connecting I = J + | with I = J —\. 
So the scattering matrix S Jl = exp(2iSji) must be a 
diagonal matrix. It is then verified that Eq. (54) reduces 
to the conventional single-channel Liischer's formula in a 
moving frame discussed in subsection II B. 

B. Two-channel case 



In this subsection, we generalize the results of a single 
channel case to the two-channel case. The formalism is 
the same as in the single channel case except that we need 
two interpolating operators (Jj(x) with j = 1,2 denoting 
two different channels. In this shown in Fig. 1, 

the two-point function has the following form, 



C(P)= / a r:q [Z 1 A 1 Z2A 2 }jg- q crl. q + / a r:q [Z 1 A 1 Z2A2}jg- q iKg h;q:ql [Z 1 A 1 Z2A2]hr-,q"jt. q ' 

J L;q J L;q,q' 



(55) 



I 

Here indices g, h, and r also refer to the channel and with the definitions 
take the value 1 or 2. In Eq.(55), we have also utilized 
the following shorthand notations: Z rJ A rJ (q) = I d i xe iq ' x ((t> r j{x)4> rJ {Q)) , (57) 



[ZiAiZ 2 A 2 



J 9-1 



<),, /,, A,, U,) Z !i:1 A !i:1 iP <!)] 



(56) 



A i; i(g) 



ha) 

■f is 



(58) 
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A*2(g) 



<7 2 



is 



(59) 



So we have two indices for the propagators and their 
dressing functions: Aj.j (and the corresponding mass val- 
ues rrij j which enters the propagator) and Zj-j. The first 
index j = 1,2 designates two different scattering chan- 
nels. The second index I = 1,2 now denotes particle 
types: I = 1 for the baryon and I = 2 for the meson. 

The kernel iig/j is again related to the Bethe-Salpeter 
kernel BS g h which now becomes a matrix in channel 
space. The matrix elements of interpolating fields Oj be- 
come vectors in channel space. We have also assumed 
that Oj carry definite parity quantum numbers. So simi- 
lar to Eq. (46), we have, 



^•(k*) = (0|^(0)|k*,-k*;i;in)| k ,| =g . 



\4(k*) = (k*,-k*;j;out|o-t(0)|0>| k .| =g; ' 
In angular momentum basis, expansion (47) becomes 



(60) 



A j {k*)=V^A j .j Ml Y l J j}(k*) 
A' j (k*) = V^A' j . JMl Y l J i(i*) 



(61) 



The factors F and M entering the quantization condi- 
tion (53) in the previous subsection have both become 
2x2 matrices in channel space. We may expand them in 
terms of spin spherical harmonics, with the coefficients 
F being diagonal in channel space: 



Fij;JMl;J'M'V = ^ijFi;JMl;J' M'V ■ 

The diagonal element is given by 



(62) 



i;JMl;J'M'V 



8tvE' 



■(Sjj'Smm'Sw + iF t 



(FV) 

i-JMl-J'M'V i 



(63) 



The scattering matrix M, however, is not diagonal in 
channel space. It is still related to the non-relativistic 
scattering amplitude via, 



M. 



ij;JMl;J' M'V 



8irE*M, 



(NR) 

ij-JMl-J'M'V 



(64) 



where M^J^ri-j'M'V * s defined in Eq. (A18). Accord- 
ing to FIG. 11 (d), Eq. (45) still holds except that the 
quantities involved have all become matrices or vectors 
in channel space. The poles in C FV still yield the de- 
sired quantization condition (53) with the understanding 
that both M and F have now become 2x2 matrices in 
channel space. 

We are now in a position to write out the quantiza- 
tion condition to a form that is comparable to the con- 
ventional multi-channel Liischcr formula obtained in the 
previous section. For this purpose, we define a matrix 



U 



(FV) 



using the matrix in Eq. (63) 



U 



(FV) 



'FF V) +i 



i;JMl;J'M'V 



F 



(FV) 



(65) 

JMl;J'M'V 

Now consider an arbitrary moving frame. As we will 
show in appendix B, sec Eq. (B21), the matrix element 



F 



(FV) 



i;JMl;J'M'l 



IV is in fact related to the corresponding ma- 



trix element M^-jmi-tm'V that have appeared in our dis- 
cussion in subsection II B, c.f. Eq. (29): 



r i;JMl;J'M'V 



i l 1 Mf : jMhJ'M'l' 



(66) 



This expression is valid up to terms that are exponen- 
tially suppressed in the large volume limit. Note that 
if the Hamiltonian conserves parity such that scattering 
only occurs for 1 = 1', the two matrices become identical. 
Substituting the definitions of Mij-jMi;j' M'V (64) and 
Fij:jMi;j'M'i'{62) to the quantization condition Eq. (53), 
we arrive at the following result: 



\ r l;JMl;J'M'l' " r V J ll\ r l;JMl;J' M'V 
~ S Jl( F (FV) _.x 
q' J 12\ r l;JMl;J'M'V °) 



V Y q » °21\ r 2;JMl;J'M'V 

(FV) 

2-JMI-.J' M'V 



i) 



(A 



S 22( F : 



(FV) 

2-JMl-J'M'V 



i) 



= 



(67) 



Comparing (18) with (67), we find that the forms of 
these formulae are completely equivalent although we 
have used different methods to obtain them. This con- 
clusion is valid up to terms which vanish exponentially 
with the box size. 



IV. DISCUSSIONS AND CONCLUSIONS 

In this paper, we have generalized Liischcr's formula to 
the case of multi-channel two-particle (one spinlcss, one 
spin-1/2) scattering in a cubic box. The generalization 



was done using both non-relativistic quantum mechanics 
and quantum field theory. We verified that, up to terms 
that are exponentially suppressed in the large volume 
limit, both methods yield compatible results. 

Although wc only consider the case of mcson-baryon 
scattering in this paper, using similar techniques, it 
should not be too difficult to generalize the results to 
the case of scattering between hadrons with other spin 
configurations. In fact, using similar notations as in 
Refs. [33, 34], one should be able to obtain corresponding 
formulae suitable for nuclcon-nuclcon multi-channel scat- 
tering. Another interesting direction is the correspond- 
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ing formulae in a box with different boundary conditions, 
which turns out to be useful for practical reasons. In par- 
ticular, the formulae obtained in this paper can readily 
be generalized for anti-periodic or twisted boundary con- 
ditions. 

A unique feature of the multi-channel scattering 
Liischcr's formula which differs from that in the single- 
channel case is that, the corresponding equation is not 
a one-to-one relation between the energy and the cor- 
responding scattering parameters. Therefore, even if 
one can construct appropriate correlation functions to 
obtain the two-particle energy E in lattice simulations, 
Liischcr's formulae only sets up constraints among the 
energy E and the 5*-matrix parameters S-'j(E). Further 
physical inputs are needed to really pin down these scat- 
tering parameters in a multi-channel scenario. 

Finally, let us discuss some possible applications of 
Liischcr's formula in multi-channel scattering. One typi- 
cal example is the antikaon-nucleon scattering. The scat- 
tering amplitude of antikaon-nucleon is of fundamental 
importance in the study of A(1405) resonance which just 
exceeds the scattering threshold. There is a strong cou- 
pling between A'N and T,tt channels when the energy 
exceeds KISS threshold. It becomes a problem for two- 
particle scattering, one with spin i and one with spin 
in two channels. In Rcf. [48], the authors used two chan- 
nel Lippmann-Schwinger equation to study the problem. 
In Ref. [49] , the same problem was addressed using unita- 
rized chiral perturbation theory. In principle, this prob- 
lem can also be studied using Liischer's formula in lattice 
QCD simulations, although that requires more data than 
we currently can acquire. Another example is the DN 
and 7r£ c coupled channels scattering studied in Ref. [50] . 

To summarize, in this paper we generalize Liischcr's 
formula to the case of particles with spin, beyond the 
inelastic threshold in COM frame and MF respectively. 
Using a quantum mechanical model and quantum field 
theory, a relation between the energy of the two-particle 
system and scattering matrix elements is found. It is 
verified that the two methods yields the same result if 
we neglect terms that are exponentially suppressed in 



the large volume limit. Although we focus on the case of 
scattering between two particles: one with spin and the 
other with spin \, we do not see any essential difficulties 
to generalize the situation to any spin with any number 
of channels. We hope that these relations will help us to 
study multi-channel hadron scattering between particles 
with spin using lattice QCD in the future. 
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Appendix A: Single-channel spin-dependent 
scattering 

In this appendix, we first briefly review single channel 
scattering of a particle with spin in quantum mechanics 
in infinite volume. After that, single channel Liischer's 
formulae are collected which can be found in Refs. [21, 
24, 33, 34]. Similar formulae in the multi-channel case 
arc listed afterwards at the end of this appendix. 

To avoid inessential complications, wc shall only dis- 
cuss the case of two stable particles with spin and spin 
\ in COM frame. In non-relativistic quantum mechan- 
ics, after factoring out the center of mass motions, the 
asymptotic form of the wave function is: 

^(r) = X Ie 4k ' r + J2 xl,M { s N s R} (k ■ f)^ . (Al) 

s' 

This form has the property that, in the remote past, it 
reduces to an incident plane wave with prescribed quan- 
tum numbers (linear momentum and spin) . In Rcf. [42] , 
the scattering amplitude between particles with spin 
and spin | is given by, 



4tt 



j 



all £< s '" 

;=o j=i-i m=-j 



(A2) 



where the auxiliary spin spherical harmonic function 
?)jMis(f) in spin-space is defined as 



-/.it 



2W(f) = r'xJ' - Y Mi) 



(A3) 



The dot here indicates an inner product in spin space. 
This means that the scattering matrix is diagonal in 



angular-momentum basis: 



M 



(NR) _ M (NR) 

JMl;J'M'V ~ lvx J'M'l'\JMl 

= Sj>jSm'mSi'i-—(S ji 
zik 



1). (A4) 



For single channel spin-dependent scattering, the S- 



matrix is parameterized by 



Since J can 



take two possible values, J = I ± 1/2, wc will also conve- 
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nicntly denote them as g^-'i 1 / 2 ,' = gl± = cxp(2i<5;±). cidcnt momentum k, the scattering amplitude depends 
If we have chosen the z-axis to coincide with the in- only on 9: 



M ^ R) ^=2ikH E V4tt(2/ + l)(S Jl - l)S% ;ms ,S l f M . i0s Y lm (v) 



(A5) 



;=o j=i-i 



Due to the Clebsch-Gordan coefficients Sf M . 0s and 



S JM;ms>> We find: M 



s and m = M — s' . Thus, in 
the previous equation we ignore the sum over M and m 
for a given pair of s and s' . It is then clear that M{&) as a 
matrix in spin space can also be written in the following 
form 



M(6) = f(6)+ig(e)(<T-e) 



(A6) 



where f(6) and g(9) are known as no-flip and spin-flip 
amplitudes, respectively, er is the Pauli matrices in spin 
space and e is a unit vector perpendicular to the scat- 
tering plane. With this convention, these functions are 
given by 



f(0)=J2[{l + l)fi + + lfi-}Pi(cosl 

1=0 
oo 

g(fl) = Eift i -fi ]p/Ng) — 



(A7) 



I 

with fi±(k) = e lSl± sm(5i±)/k. These arc identical to 
what have appeared in the early literature, see for exam- 
ple, Rcf. [51] and Ref. [52]. 

Let us now enclose the system we considered in a 
large cubic box with periodic boundary conditions ap- 
plied in all three spatial directions. The interaction is 
only present in the inner region while in the outer re- 
gion, the radial wave function becomes a superposition 
of ni(kr) and ji(kr), 

V(r) = X>-n(*)jH*r) + Pji(k)m(kr))Y^(r) (A8) 

J Ml 

which leads to the identification tan(<S,/;(fc)) = 
fi,n{k) / Oiji{k). The Green's functions, which are singu- 
lar periodic solutions of the Helmholtz equation, has a 
similar expansion. In this case, it is similar to those in- 
troduced in Ref. [4] (see also e.g. Ref. [33, 34]): 



Gjmi(j, k 2 ) 



± ~4 L ~ [Yjh(*)ni(kr) + ]T Mjui^MU'^Yp^j^kr) 

71 \ J' API' i 



(A9) 



Comparing Eq. (A8) with Eq. (A9) and following similar 
steps as in the derivation of the conventional Liischcr's 
formula, one finally finds 

det[tanSji(k)M jA4i-r An' - 5jj'5mm'8w\ = ( Al °) 



I 

which is the same as in Ref. [33]. The explicit form of 
■MjMl;J'M'l' is given in Ref. [33] which we quote here: 



Mjmi;J'm>i> = E M lm ., rm >(lm\s\JM){l'm')-s\J'M') , 



(All) 



M 



Imd'm' 



(-i) 



i w 



E Etttt^ 1 ;^ 



m : tn,l' m' 



(A12) 



with k = ki/(27r) and the zeta function Z tn (\;n 2 ) and the coefficients Ci m ^ n ,i'm' ar e given respectively by 



J 



z tn {\-K 2 )= 

n£% 3 



n 2 — k 2 



(A13) 
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C lm , tn ,v m > = i l ~ t+l 'J ^ 21 ^ (I0t0\l'0)(lmtn\l'm') 



(A14) 



r 



The corresponding formulae in moving frames have 
been obtained in Ref. [21, 24]. For example, instead of 
Eq. (All), we have, 



-1)' 



I m.tn.l' m' 



(A15) 



t ,71 



with k = kL/(27r), the parameter 7 is the Lorentz boost 
factor associated with the moving-frame and the zeta 
function Zf n (l\K 2 ) is defined by 



Z? n (l;K 2 )= 

nePd 



|n|*y t „(n) 



(A16) 



the coefficients Ci m .t n .i'm' are the same as in non-moving 
frames, e.g. Eq. (A13). In the above formulae, 



V A = \v 



r = 7 : (n+ -ad) nef 



(A17) 



where a and d are given as in subsection II B. Then 
Eq. (All) still holds except that one has to use the 
moving-frame versions (M j M i;j>Mn> and Mf m . Vm ,) in 
place of the non-moving frame versions. 

Extension of the above formulae to the case of multi- 
channel case is straightforward. For example, Eq. (A4) 
is generalized to 



M 



(NR) 

ij-JMl-J' M'l' 



- M (NR) 

~ 1U ij;J'M'l';JMl 

= SjijSm'mSp 



2z ^/ kj^k 



(A18) 



Similarly, Eq. (All) is modified to: 
Mi,j Ml ,j, Mn ,= M% m . Vm ,{lm l -s\JM){l'm l ^s\J l M' 

mm' s 

(A19) 

with Mf. lm . lljn , given in terms of the zeta function 
Z$ tn (l; Hi 2 ) like Eq. (A15), with m = kii/(27r). 



where the function /*(k*) has no singularities for real 
k* and falls off fast enough at k* — >• 00 so as to ren- 
der the summation convergent. To simplify the matter, 
we assume that it is spin- independent and thus can be 
expanded into spherical harmonics as in Ref. [45]: 



(B2) 



One would like to study the behavior of S(q*) in the 
limit of large L. For a fixed L, if there is no term in 
the sum with k* 2 ~ q* 2 , we can replace the sum by an 



integration. The singularity at k* 



q* 2 for large L 



forbids this simple replacement. Basically S(q*) will split 
into two parts, one of which can be approximated by 
the principle-valued integral, the other being the finite 
volume correction. This has been established in Ref. [45] 
and we directly quote their final results: 



S(q*) 



d 3 k* /*(k*) 
(2tt) 3 q *2 _ fc *2 



^2fhnC[ m {q* 2 ) , 



(B3) 



E 



V -k* l V^Y lm (k*) 



L 3 iDk q* 2 — k* 2 



-V 



d 3 k* e ai -^~ k '^ 
(2tt) 3 q* 2 -k* 2 ' 



(B4) 



where V stands for the principal-value prescription. This 
summation formula will be utilized shortly. 




Appendix B: Calculation the the kinematic factor of 
loop integration/summation in a single channel 

In this appendix, we use the notation: k* = |k*|, q* = 
|q*| unless otherwise stated. 

The generic finite-volume corrections that we are in- 
terested in has the following form: 



w=Ae < r(k * } 



L 3 ^ Wk q* 2 — k* 2 



(Bl) 



FIG. 2. The dash line represents the meson propagator, and 
the solid line represents the baryon propagator. A circle with 
a symbol iK inside stands for the Bethe-Salpeter kernel which 
consists of all amputated two-particle irreducible diagrams. A 
circle with a symbol a or denotes the interpolating opera- 
tors in Eq. 33. 

We now come to the correlation function C (P) defined 
in Eq(33). It can be expressed in terms of the Bethe- 
Salpeter kernel iK through the series shown in Fig. 2. 
The loop integration/summation appearing in the figure 
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has the following form: 



k 



dko 



/(fc 0) k)(fc"7 M + mi ) 



2tt (fc 2 
where fc = (feo,k), P = 



m 2 + is)((P — fc) 2 — to 2 + ie) 

(B5) 

(E, P) being the corresponding 
I 



four-momenta and mi and m 2 being the masses of the 
baryon and the meson respectively. The function /(fc) 
contains the energy-momentum dependence arising from 
the kernels as well as that from the dressed propagators. 
The properties of /(fc) is such that there exists no singu- 
larities for real k and its ultraviolet behavior is to render 
the integration or summation convergent. Integrating 
out fco one gets 



I = — 



L 3 

k 



/K k ,k)c(k) 



2w lk ((E 



wik) 2 - w 2 



f(E + w 2kl k)D(k) 



2w 2 k((£ + W2k) 5 



w 



Ik' 



(B6) 



If we assume /(k) is an even function for k 

f C(k) = (fc 7o+TOi)| fco=t „ lk 
[ D(k) = (fc° 7o + m x ) \ ko =E+w 2k 



I 

contain the two-particle finite- volume poles is Ji, 
_/K,k)C(k) 



2w lk ((E - wi k ) 2 - w 2 k ) 



(B8) 



(B7) 



and w\k = yk 2 



It is seen that the two-particle pole singularity in I\ is 
located at E = + W2k- To determine the finite vol- 
ume correction in more detail, we express the term 1\ in 
another form by transforming it in to the COM frame. 



W2k = \/(P k) 2 + m 2 are the In COM frame the two energies are: 



w 



two energies. Note that both C(k) and D(k) and thus 
the integral / are matrices in Dirac space. However, as 
mentioned in the main text, normally the interpolating 
operators that we formed to create the meson and baryon 
carry definite parity. This means that, for practical ap- 
plications, 70 takes its eigenvalues, ±1 times a unit ma- 
trix (in spin space). Thus the integral I is written into 
two parts 1\ and I2 corresponding to the two terms in 
Eq(B6). The second term I2 does not contain the finite- 
volume singularities in the kinematic region of interest 
and therefore can be replaced by the corresponding in- 
tegral in the large volume limit. The term which does 



lk 



*2 



w 2k 

tain: 



h = - 



as given in section II B. Thus we ob- 



f(k*)C(k*) 



L 3 ^ 2w lk ((E* - W * k ) 2 - w 



2 ) 

2k' 



— Y 

2L 3 ^ 



Kk /*(k*)C(k*) 
Wik q* 2 — fc* 2 



(E* 



u 2k 



AE* 



■w 



lk 



(B9) 



By using the summation formula Eq. (B3) we men- 
tioned at the beginning of this appendix, we obtain: 



h = {-i)V 



d 3 k* /*(k*)C(k*) 
J2^f q* 2 -k* 2 



(E*+w* lk ) 2 - W * k 



8E* 2 w* lk 



2E* 



111) 



(BIO) 



I 

In order to write I\ as the infinite-volume result together tion in the above formula by a Fcynman +is prcscrip- 
with a correction, we replace the principle- value integra- tion in the propagator and a "delta-function" term which 



J 



picks out the I = part of /*. 



h 



1 f d 3 k* /*(k*)C(k*) (E*+w* lk ) 2 ~w* 2 2 , q*fSo(q*)C(q*) *C(q*) 



2E* J (2tt) 3 q* 2 - fc* 2 + ie AE*w 



lk 



8ttE* 



2E* 



Finally, we arrive at our final result for /, 

I = loc + IFV , 



(B12) 



r 



where 1^ is the infinite volume result for the original 
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loop integral with the appropriate Fcynman's prescrip- 
tion, = I% t00 + li- As mentioned earlier, 1^ contains 
no finite volume singularities. The two-particle singular- 
ities are contained in the finite volume correction term, 
Ifv which is given by 

lm 

(B13) 



where V stands for the principal-value prescription. 

In Rcf. [47], for the scattering of two particles 
with different masses, a relation between Cf m (q* 2 ) and 
Zf m (l;K 2 ) is found, 

with k = q*L/(2n), Zf m (l;K?) is defined by (A16). 
This equality holds up to terms which vanish exponen- 
tially with the box size. Thus, according to the expres- 
sions for F$£ l .j, MIV and M^ m .j, M , v i.e. Eq. (B16) and 
Eq. (A19) we have the following equality: 

F JMl;J'M'l' = ■M-JMl;J'M'l' ( B19 ) 

which is also valid up to terms that vanish exponentially 
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